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Abstract 



We consider interference effects between p-wave resonance scattering ampli- 
tude and background s-wave amplitude in low-energy neutron scattering from 
a heavy nucleus which goes through the compound nucleus stage. The first 
effect is in the difference between the forward and backward scattering cross 
sections (the Pi • P/ correlation). Because of the chaotic nature of the com- 
pound states, this effect is a random variable with zero mean. However, a 
statistical consideration shows that the probability distribution of this effect 
does not obey the standard central limit theorem. That is, the probability 
density for the effect averaged over n resonances does not become a Gaus- 
sian distribution with the variance decreasing as n~ 1 / 2 ("violation" of the 
theorem!). We derive the probability distribution of the effect and the limit 
distribution of the average. It is found that the width of this distribution does 
not decrease with the increase of n, i.e., fluctuations are not suppressed by 
averaging. Furthermore, we consider the cr- (pj x pA correlation and find that 
this effect, although much smaller, shows fluctuations which actually increase 
upon averaging over many measurements. This behaviour holds for e > T p 
where e is the distance to the resonance, and T p is the resonance width. Lim- 
its of the effects due to finite resonance widths are also considered. In the 
appendix we present a simple derivation of the limit theorem for the average 
of random variables with infinite variances. 



PACS: 24.60.-k, 25.40.Dn, 05.40.-a, 02.50.-r 
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I. INTRODUCTION 



Due to chaotic nature of compound nuclei, positions of s and p resonances in neutron 
scattering from a heavy nucleus, and amplitudes involving these states, are uncorrelated. 
This gives rise to an unusual statistical effect in the asymmetry of the transmission of 
neutrons with positive and negative helicities [0J . This asymmetry corresponds to the a ■ p 
correlation. It violates parity conservation, and is produced by the weak interaction in the 
nucleus, which mixes the s and p neutron partial waves. The magnitude of the asymmetry 
is strongly enhanced if the neutron energy is tuned into the p resonance f2j. In this case its 
magnitude is determined by the perturbative mixing 



of the s and p resonances by the weak interaction W. The matrix element between the 
compound states behaves as a Gaussian random variable, and 77 is also a random variable 
with zero mean. The characteristic mixing (and asymmetry) can be estimated simply as 
T]c ~ w/D, where w is the root-mean-square matrix element, and D is the mean spacing 
between the compound nucleus resonances. However, when one takes a closer look at the 
probability distribution of 77, it turns out that its variance is infinite! This effect originates 
from a high probability to find small spacings E s — E p , which results in the slow decrease of 
the probability density, f(rj) oc I/77 2 . 

The standard Central Limit Theorem (CLT) is not applicable to such random variables. 
Instead, it can be shown that if we consider the statistics of the average of n such variables, 
the probability density of the sum (rji + . . . + r] n )/n becomes independent of n at large n, 
i.e., fluctuations of rj are not averaged out ||. This contrasts the usual situation where CLT 
would give a Gaussian distribution whose width decreases as 1/ y/n. This unusual behaviour 
is explained by the fact that among n uncorrelated rji there is a large probability to find 
one, whose magnitude is n times greater than their typical value, rji ~ nr] c . Such rji will 
always dominate the sum and ensure that fluctuations are not suppressed by averaging. 
Another physical instance in which rare events dominate the distribution is seen in "Levy 
flights" in the random force diffusion model |3J. In a usual homogeneous system, diffusion 
is modelled by Brownian motion, where the distribution of displacements is Gaussian. But 
in the random force diffusion model, disorder induces rare but large displacements which 
dominate the distribution (the Levy distribution). Under certain values of parameters, these 
have statistics similar to those found in nuclear scattering problems. 

In this paper we analyse other effects in neutron scattering that have such unusual 
statistics. They are more conventional than those discussed in Ref. |IJ, since they do not 
involve the weak interaction. For example, in Sec. |TJ we consider the difference between 
the differential cross-section in the forward and backward scattering, due to the interference 
of the p-wave resonance scattering and the background scattering amplitude, for a spinless 
particle. This effect can be described as the Pi • p/ correlation, where p» and p/ are the 
momentum of the incident and emitted particles, respectively. We derive the statistics of the 
observable and the way it behaves upon averaging over many different nuclei, in a similar 
fashion to what was done in |IJ. We also discuss the limit of this effect when finite widths 
of the resonances are taken into account. In Section [ITT] we rederive the p« • P/ correlation 




(1) 



7/ 



E s — E p 
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for particles with spin (e.g., neutrons) incident upon a spinless nucleus, as well as study 
the effect of a different correlation between the neutron spin er and the scattering plane, 
a ■ (pj xpf). This is shown to have different statistics to the first correlation, and we derive 
a limit theorem for the average of the second correlation effect (details can be found in the 
Appendix). As it turns out, fluctuations of this average increase upon averaging, because 
when larger sets of data are considered there is a finite probability of finding an effect whose 
magnitude is ~ n 2 times larger than its typical value. 



II. CROSS SECTION ASYMMETRY FOR A SPINLESS PARTICLE 

Let us first study the simple case of the p« ■ P/ correlation for scattering of a spinless 
particle. Here we consider the difference between forward and backward elastic scattering 
differential cross sections near threshold, due to interference of the p-wave resonance scatter- 
ing with the background s-wave amplitude. For a spinless particle the scattering amplitude 
at low momenta is written using the Breit-Wigner formula as (see, e.g., ||) 

a r (n) 

f{9) = -A - % t* — costf, (2) 

JK ) 2k E - E p + \ r p ' K J 

where A is the s-wave scattering length, k and E are the wave number and energy of the 
projectile, and g p , E p , and T p , are the statistical weight, energy, capture (or elastic) 
width and total width of the p resonance, respectively. We assume that at energy E the 
s-wave background is nonresonant, and there is a p-wave resonance nearby, a condition 
which would favour larger asymmetries. This leads to expressions for the forward (+) and 
backward (— ) scattering amplitudes 

-p(n) 



where e = E — E p is the distance to the p-wave resonance, and we assume that it is greater 
then the resonance width, i.e., e ^> T p . 

The relevant observable is the asymmetry 



_ (da/dn)+ - (da/dtt)_ 
X ~ (da/dtt) + + {da/dti)^ ^> 

where (da/dQ)± are the forward and backward scattering cross sections. Substituting am- 
plitudes (^|) and taking into account that at low momenta the contribution of the p wave is 
much smaller than that of the s wave, we obtain 



r (n) 

(n) 



where (3 = Ak/g p . A typical value of this asymmetry is g p T p /AkD 
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A. Statistical analysis 



We would now like to obtain the probability density for the observable x. The capture 
width is proportional to the square of the capture amplitude. The capture amplitudes for 
complex compound nuclear states have a Gaussian distribution |J, hence, the widths 
are distributed according to the Porter-Thomas law 

where 7 = for convenience and 7 is the mean width 7 = J 7^(7)^7. 

For a given energy E the distance to the nearest p resonance in a compound nucleus is 
random. If the relative positions of the p resonances were uncorrelated it would be described 
by a Poissonian distribution 

f D (e)=D- 1 ex P (-^pj (7) 

where D is the mean spacing between the p resonances and \e\ = D/2. Correlations between 
the positions of compound states of the same symmetry, often referred to as level repulsion, 
modify the above distribution. These correlations are described by the random matrix theory 
0, and can be approximated by the Wigner law. In this case the distance to the nearest p 
resonance has the following probability density [T|]: 

.-1 ' 



f D (e) = D- 1 exp^-—j . (8) 

To avoid confusion we should stress that here e is the distance to the nearest p-wave resonance 
and not the interval between the p-wave resonances. Therefore, Eq. (H) differs from the 
Wigner-Dyson distribution. The difference between Eqs. (|7|) and (||) is not very important 
for our consideration, as long as fo( £ ) remains finite (and equal to D^ 1 ) at e — > IQ. 
Using Eqs. (|5|) and we calculate the distribution of the observable x as 



/(*) 



f o d ^f def D (e)g(-f)S (^x - -^J 



^/2nxo\x\ Jo V 2x 



v^exp ( -— - Tit 2 ) dt, (9) 



where xq = \j/(3D\ represents some typical value of the effect. The integral in @ can be 
given explicitly in terms of the parabolic cylinder functions D v [z), 

_ (27r)- 3 / 4 r(j) / X 2 \ / \X\ \ 

{X) ~ v^^N eXP \327ra;o; ^ UoV^F/ ' 

On the other hand, one can easily find the asymptotic behaviour of the probability density 
at \x\ ^> x directly from Eq. @: 

f{x) ~ x /x 2 . (11) 

The CLT will not apply to distributions with this asymptotic form, since they do not have 
a finite variance: J x 2 f(x)dx — > 00. 
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B. Limit theorem for the first correlation 



Suppose the forward-backward asymmetry is measured in an experiment where a number 
of different nuclei with similar-sized cross sections are involved. The measurement will 
then yield some average asymmetry, and we want to find the probability distribution of 
it. Otherwise, one may just analyse the asymmetries measured separately for a number of 
nuclei. Let us then consider the average of n independent random variables Xj introduced 
above, 



X = 

n 



1 n 



x=l 



In Ref. [IIJ a derivation of the limit theorem for distributions with asymptotic behaviour 
(PI) was presented. A general solution of the problem of limit distributions of sums of 
independent variables with an infinite variance for which f(x) oc can be found in 

Ref. (a > to keep the total probability J f(x)dx finite). A simple derivation of the 
limit theorem for such distributions is given in the Appendix. 

The random variable x has a symmetric probability distribution, f(—x) = f{x). In this 
case for a = 1 the limit distribution is obtained from Eqs. ( |A20| ) and ( |A22|) with a = and 
c = ttxq. So, for n — > oo the probability density F n (X) approaches its limit form 

F ^ = Iwhrr (13) 

where X c = ttxq [[[]. This is called the Cauchy distribution, and its main property is that 
it is independent of n, in particular, it does not become narrower as n increases. Therefore, 
fluctuations are not suppressed by averaging. Compare this with the standard central limit 
theorem, where the width decreases as o n = a\j y/n. 

The parameter X c for our physical observable is given by X c = ix^ / f3D = ng p Tp^ /AkD. 
Throughout the derivation we considered the scattering length A as a constant. Indeed, 
if the energy E does not coincide with an s-wave resonance, the s-wave amplitude —A 
represents the potential scattering amplitude. It does not vary strongly between isotopes, or 
nuclei of similar masses, because the nuclear potential does not vary much. The energy scale 
of its variation is MeV, similar to single-particle energy level spacing. Contrast this with the 
scale of the compound resonance spacings which are of the order of 10 eV. This difference in 
energy scales means that the compound resonance states can be treated statistically, while 
the scattering length is treated as a constant. 



C. Influence of the resonance widths on the statistics of the average asymmetry 

The above calculations have been based on the assumption that e ^> T p , so that the 
possible effects of the resonance widths have been neglected. This is justified, as long as the 
probability of finding a very small interval e ~ T p is indeed small. However, it is easy to see 
the role of the width as we increase n. As we explained in the Introduction, the nonvanishing 
fluctuations of the average depends on having one value of the effect large, Xi ~ nx c , where 
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x c is a typical value. Indeed, we can expect that if we make n measurements then at least 
one will have an energy spacing of the order e ~ D/n, thus giving a large asymmetry (^j). 
The energy denominator, however, can not be made arbitrarily small, and as smaller E — E p 
are considered, it will reach a limit \e + iY p /2\ ~ T p . Thus D/n ~ T p determines the 
largest values of n, beyond which the Cauchy distribution of the average effect begins to 
turn into a Gaussian one. Hence our statistical analysis is valid until n ~ D/T p (in heavy 
non-fissionable nuclei D/T p ~ 500). 

If we continue to take measurements after this and further increase the number of mea- 
surements n, the maximum value of the asymmetries will stabilise, being of the order 
Xi ~ XqD/T p = Tc^/PTp. Hence we will no longer have increasingly large values of the 
effect to continue the "non- vanishing" averaging. Thus, when n ^> D/T p the Gaussian 
statistics take over, the standard central limit theorem applies, and the usual l/\/n sup- 
pression of fluctuations takes place. Note that xqD/Y p 3> x in fact determines the true 
finite, but large, variance of x. Beyond this value f(x) decreases faster than x~ 2 , and it 
effectively determines the lower and upper limits in the variance integral J x 2 f(x)dx. 



III. SPIN ONE-HALF PARTICLE CORRELATIONS 



Let us now consider scattering of low-energy spin-^ particles, looking at both the p» • p/ 
and cr- (pj x pj) correlations. Again we assume that there is no nearby compound nucleus s- 
wave resonance and that the asymmetry is dominated by one nearby p-wave resonance. This 
is justified because the statistics of the average effect X for large n, F n (X), is determined 
by large values of the individual effects, i.e., by the asymptotic large- x behaviour of the 
probability density f(x), cf. Eq. (0) §. 

Consider scattering of the neutron, s = \ from a nucleus with spin /. The total angular 
momentum of the p-wave neutron is j = 1 + s and the total angular momentum of the 
compound resonance is J = I+j. The amplitude of p- wave resonant scattering at arbitrary 
angle 6 between the incoming and outgoing particle can be written as (see Ref. 0) 



x E - Ep + i Tp C ^ C t^ v^Un.) Vll? (E), (14) 

where I z , n&, s z , and I' z , n' k , s' z , describe the projection of the target spin, the direction of 
the neutron momentum and the projection of the neutron spin in the initial and final states, 
respectively, T p n ^ is the capture width for the neutron with angular momentum j, the Y\ m 
are the angular wave functions, and C T f z -- are the Clebsch-Gordon coefficients. 

Consider scattering of a neutron incident along the x direction off a spinless (J = 0) 
nucleus. We quantise the neutron spin in the z direction and consider neutrons scattered 
in the xy plane. The s-wave scattering amplitude is simply — A5 SzS ' z . Having in mind that 
we need to calculate interference terms between the s and p waves, we can write the p-wave 
scattering amplitude ( |14"D in the following form 

fM = ~Yk E \ C t^J Y M) E _ ^ +iv Y lm {n k )^{E). (15) 



6 



A. First correlation 



The resonance p-wave amplitude for forward and backward scattering (see Sec. |TJ) with 
= § is 

f± = T — ^ (16) 

which is similar to the spinless particle scattering (Eq. ^), with the parameter g p = 1. 

3 
2 



Similarly for j — ~ states the amplitude is 



/± = T — P3/2 . (17) 

Jp3 / 2 2k E - E p + ±r p V ; 

which is similar to spinless particle scattering with g p = 2. This means that the statistics 
derived for the spinless particle Pi • p/ correlation in Sec. [H] are valid for the case where spin 
is included. In fact, since we do not know whether the nearest resonance is pu2 or p 3 / 2 we 
must combine the two distributions. 



B. Second correlation 

The second correlation a ■ (pj x pj) between the direction of the spin relative to the 
scattering plane is, of course, specific to particles with a non-zero spin. To calculate the 
asymmetry of the cross section with respect to flipping the spin, we take the initial neutron 
momentum along the ^-direction as before, and look at the difference between the scattering 
amplitude in the +y direction, / + , and that in the — y direction, /~. Equation (|15D yields 
the p-resonance scattering amplitudes in the +y and — y directions for j = | 

. p(n) 

h ^ 2kE-E p + ^T; 1 ; 



and similarly for j > = ~ we obtain 



i 



r (n) 
1 P3/2 



Jp ^ ^2kE - E p + iT p V ' 

which differs from the j = \ case only by sign. 

Thus the total scattering amplitude for the cr ■ (p f x p^) correlation is given by 

p(n) 

f± = _ AT v UL ip . (20) 

J 2k E — E p + |T p 1 ; 

where = — 1 for j = | and r) p — +1 for j = |. Then, taking into account that the second 
term in Eq. (0), which represents the p-wave contribution, is much smaller than the first 
one, we obtain for the observable difference of the corresponding cross sections [see Eq. (H)] 
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As we discussed above, the scattering length varies weakly. The same is true for the total 
width of the compound resonances T p . Its fluctuations are small because it is dominated by 
the radiative width, given by the a sum of a large number of partial widths due to transitions 
into all lower- lying nuclear states. Introducing (3 = 2kA/r] p T p , and taking into account that 
e = E — E p is usually much larger than the resonance width, |e| ^> T p , we obtain for the 
asymmetry 

r (n) 

The typical size of this effect f] p r p n ^r p /AkD 2 = (j] p T p n ^ /AkD)(T p / D), is much smaller than 
the first correlation, by a factor of T p /D. However, this observable has a e~ 2 dependence on 
the distance to the nearest p resonance, while the first correlation was proportional to e . 
The e~ 2 singularity emphasizes even stronger the possibility of small denominators. Note 
also, that for a given scattering length A the sign of this interference effect is always the 
same. We will see that this leads to a very different statistics of the a ■ (p« x pf) correlation. 



1. Statistical analysis 

Let us derive a probability distribution for the observable x given by Eq. (p2|). Similarly 
to Sec. [II A| we have 

where we again use 7 for the capture width Tp , and the probability densities #(7) and 
/d({t) are taken from Eqs. (|]) and (f|), respectively. Assuming that the scattering length is 
positive, hence, (3 > 0, we calculate the above integral and obtain 



^TTX(X + TCXo) 

and f(x) = for x < 0, where 



f(x)= * X ° -, x>0, (24) 



" = W (25) 



characterises typical values of the asymmetry ( j22[ ). The asymptotic behaviour of this prob- 
ability density at x ^> xq is 



The probability density f(x) is normalized as J °° f(x)dx = 1. However, the corresponding 
mean value J f(x)xdx is infinite, and the integral for the variance J f(x)x 2 dx diverges even 
faster than that for the first correlation [cf. Eq. (fLl]) ]. This signifies even larger fluctuations 
of the second correlation effect. 



S 



2. Limit theorem for the second correlation 



The probability distribution of the second correlation (p6| ) corresponds to a = ~ (see 
Appendix). Using the asymptotic parameters C\ = and c 2 = \/x /ir [compare Eqs. fl2"o| ) 
and ( |A1|) 1 we obtain c = v^Co an d /5 = 1 from Eqs. ( |A9] ) and ( |A10|) . In fact, it is possible 
to calculate the Fourier transform of f(x) of Eq. (E3J) explicitly, 



f(w) = e i7TX0UJ 7r- 1/2 T(l/2,iiix Lu) 

= e inX0U} [l-v^(l± 1/2 + 0(u 3 / 2 )] , (27) 

where T(. . .) is the incomplete T-function, and ± corresponds to u ^ 0. 

It follows now from Eqs. (|A~20l) and ([A24]) that the limit distribution F n (X) of the 
average effect X is given by 

/jiT-n p -nxo/X 

F n (X) = ^ e -^, (X>0). (28) 

This equation shows explicitly that as the number of effects included in the average increases, 
the distribution widens proportionally to n. Accordingly, the typical values of the average 
also grow as X ~ nx . 

To understand this recall that the second asymmetry is inversely proportional to e 2 . 
Thus, when we consider the average of n such variables, one of them is likely to have 
e ~ D/n, which makes it n 2 times greater than the typical value xq. This variable will 
dominate the average and give X ~ hxq. Also, to describe a possible experiment more 
realistically, one must combine the distributions with different r] p remembering that the 
probability to find a close P3/2 resonance is twice that of a pi/ 2 resonance. 

The role of resonance widths can be understood in the same way as it was done in 
Sec. [11 C| . The limit statistics ( |28| ) is valid until n ~ D/T p , i.e., when the characteristic 



value of the maximal effect that dominates the average X requires denominators as small as 
e ~ D/n ~ T. When n 3> D/T p the usual statistics of the Central Limit Theorem become 
valid, and we eventually have typical values of the average decreasing as \/\fn. 



IV. CONCLUSION 

We have considered interference effects between p-wave resonance neutron scattering 
amplitude and background s-wave amplitude in compound nuclei and found that these 
effects do not obey the Standard Central Limit Theorem. That is, the probability density 
of the average effect over n measurements, X = - YH=i x i > does not tend to a Gaussian 
distribution with variance a 2 = o\/n for large n. We have examined two effects, (i) the 
Pt Pf correlation, which corresponds to the forward-backward asymmetry of the differential 
cross section, and (ii) the er • (pj x p^) correlation, which describes the asymmetry of the 
scattering cross section with respect to flipping the spin relative to the scattering plane. 

The first of these was found to have a distribution with asymptotic behaviour f(x) oc 1/x 2 
for large x. In this case the limit theorem for the distribution of the average X tends to a 
Cauchy distribution F n (X) = X c /[ti(X 2 +X 2 )]. This is independent of n, so the typical value 
of the average (X c ) does not decrease with increasing number of measurements. Physically 
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this is understood by the following arguments. The asymmetry in question is inversely 
proportional to the spacing between the incident neutron energy and the energy of the 
closest p-wave resonance (x oc s^ 1 ). If we have n measurements, we have a high probability 
that one of these spacings will be of the order e ~ D/n where D is the mean p level spacing. 
This will produce an asymmetry of the size x ~ nx where xo is a typical value of the effect. 
Thus the typical average value is X ~ xq, nonvanishing with increasing n. 

The second correlation we considered produces a much smaller effect than the first cor- 
relation. However, it has been found to have a e~ 2 dependence, giving a distribution with 
asymptotic behaviour f(x) oc x~? for large x. This means that there is a higher proba- 
bility to obtain relatively large values of x, compared to the first correlation. As a result, 
typical values for the average of n asymmetries actually increase with increasing number of 
measurements nasX~ nxo- 

Above we assumed e > T p where e is the distance to the resonance and T p is the resonance 
width. When we consider the the influence of the resonance widths, it is found that they 
affect the distribution by limiting the size of the effect, x. Indeed, the minimum value for 
the denominator \e + iT p /2\ is given by ~ T p , hence, the maximal possible effects are limited. 
We have found that our analysis of the statistics of the averages X is valid for n < D/Y p 
(for heavy non-fissionable nuclei D/Y p ~ 500). For n ^> D/T p we expect the average to 
once again obey the standard CLT and vanish with increasing n. 

Because of the interest in scattering problems, it would be of benefit to actually perform 
the experiments discussed in this paper. If one measured the first correlation in low-energy 
neutron scattering off different isotopes of heavy nuclei, then we expect to see an observable 
effect in the average. This would not decrease until the number of measurements n ^> 500. 
(Note, however, that the total number of relatively stable isotopes ~ 1000). Because it is 
hard to measure the forward and backward scattering, the experiment would have to be 
performed at some angle 6 to the axis of incidence. This would only cause an extra factor 
of cos 6? in the value of the effect. 

The second correlation is much smaller (by a factor T p /D ~ 10~ 3 ). This would make 
it much harder to observe, but it would have an increasing typical value upon averaging 
over many measurements. This means that one might be able to observe the effect after 
performing the experiment over many isotopes. 

APPENDIX A: LIMIT THEOREM FOR PROBABILITY DISTRIBUTIONS 

WITH INFINITE VARIANCES 

Consider a random variable whose probability density has the following asymptotic be- 
haviour: 




(Al) 



with < a < 2, and is normalized in the usual way, J f(x)dx = 1. The existence of 
the mean, J xf(x)dx, depends on whether a is greater or less than unity, but the variance 
integral J x 2 f(x)dx is infinite in both cases, and the standard Central Limit Theorem is 
inapplicable. 
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To derive the limit statistics of the average X = - Ylx=i x « °f n independent random 
variables x^, we use characteristic functions (or Fourier transforms) 



= / e-^f(x)dx. (A2) 

J —oo 

The Fourier transform of the probability density F n (X) of the average X is given by 

/OO f I 1 n \ n 

e~^ x dX / 5 X - - *i II 
^ V n i= i J i= i 

n 

= n/(w/n)= [/>/n)] n . (A3) 
i=i 

Thus, the form of F n (u) for large n is related to that of /(a;) at small a;. This is in turn 
decided by the large-x asymptotic behaviour of f(x) given by Eq. (|A1|) . 

For 1 < a < 2 the integral J xf(x)dx converges and /(a;) can be written as 

f{uj) = l-iu J xf{x)dx + J (e~ iuJX - 1 + iujx)f(x)dx. (A4) 

Let us consider the contribution of the interval from to +oo to the last term above. Using 
the asymptotic form ( |A1| ) we present it as 



e~ wx -l + iux) 



C2 



CO 



dx 



f( x ) -^£ri dx + c 2 I (e— - 1 + ^*)— . (A5) 



If we assume that f(x) approaches its asymptotic behaviour sufficiently rapidly, e.g., \f(x) — 
C2/x a+1 \ oc 0(l/x a+2 ), then the first integral above behaves as 0(uj 2 ) at uj — > 0. To calculate 
the second integral we turn the integration path into the complex plane by changing the 
variable uox = it (for uj > the t is real positive), which gives 



c 2 e l -u a / r a - x (e-* - 1 + it)dt, (A6) 

where the integral is a representation of the T-function on a segment of the negative argument 
axis, r(-a) = -T(-a + I) /a 0. 

Applying the same procedure to the integral over (— oo, 0) in expression ( |A4j ), and turning 



the integration path into the complex plane by using tux = —it (for to > and real positive 
t), we obtain the expansion of f{oj) at small uj: 

"X / \ / - 7TQ! \ „ ( CX 1 ) , , 

f(uj) = l-iuja - (c ie ~ l - + c 2 e 1 -) uj a — + ... A7 

a 

where a — J xf(x)dx is the mean value. The expansion for negative uj can be obtained from 
the above by simply replacing uj with and complex-conjugating the exponential phase 
factors in the second term. Finally, at the same level of accuracy, we can re-write expansion 
( |A7| ) in the form valid for positive and negative small uj: 
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c I 1 + ip sign (a; J tan — 



where 







(ci + c 2 

c 2 - Ci 
c 2 + cl 



T(l — a) TTOi 

i cos — , c > 0, 

a 2 

, -1 < /5 < 1, 



(A8) 

(A9) 
(A10) 



and sign (a;) = ±1 for u > and < 0, respectively. The parameters c and (3 [TIJ are 
determined by the asymptotic behaviour ( |A1|) of the probability density. The value of (3 
depends on the asymmetry of the probability density f{x). 

The final form ( |A8| ) is very convenient. If we consider a random variable Xi shifted 
with respect to x, X\ = x — a (a is an arbitrary number here), its characteristic function 
would differ from that of x by a simple phase factor, f\(uj) = e luja f(uj). On the other 
hand, the asymptotic behaviour of the probability density, Eq. ( |Al|) is not affected by this 
transformation. Therefore, the phase factor in Eq. ( |A8[ ) can always be eliminated by this 
simple shift. 

For < a < 1 in Eq. (|A1|) we re-write the Fourier transform as 



I +- / {e~ lux - l)f{x)dx 



The contribution of positive x to the above integral can be transformed into 

C2 



(e-^-1) f(x)- 



a+1 



dx + c 2 I (( 





00 dr 

-twx _ 



X 



a+1 



(All) 



(A12) 



Provided the difference between f(x) and c 2 /a; Q+1 decreases as x a 2 or faster, as x —>■ +oo, 
the first integral can be expanded in powers of cu, with the leading term given by 



lit) 



x 



c 2 



X 



a+1 



dx . 



(A13) 



The second integral in ( |A12| ) is transformed by variable substitution oox = it (for oo > 0) 
into [cf. Eq. flAJ)] 



cae'^w" / r a-1 (e~* - l)dt, 
'o 



(A14) 



which again gives the T-function ||12|| . After we apply the same procedure to the negative-x 
part of the integral in ( [A12j ), the expansion of f{uj) at small uj is established in exactly the 
same form as that of Eq. (|A7|) (for uj > 0). However, for < a < 1 the parameter a is no 
longer the mean value. Instead, it is given by 



x 



fix) 



Cl 



\x 



a+1 



dx + I ./• 





X 



a+1 



dx . 



(A15) 



Also, the next term omitted in Eq. ( |A7|) may now be greater than 0(u 2 ). Nevertheless, the 
small-u; behaviour of the Fourier transform is still represented by Eq. (|A8|) . 
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If a = 1 in Eq. (|Al|), the expansion of f{uj) also contains ujlnu terms. In this case it 
can be presented as 



f{oo) ~ 1 — iua — c\uj\ 



1 — i—(3 sign(cj) In \uj\ 



71 



(A16) 



where c = |(ci + C2), which can be obtained from Eq. 
( PUD , and 



) at a — > 2, /3 is given by Eq. 



(C2-Ci)(l-C)+ / 

J — c 



/(*) 



Cl 



1 + x 2 



dx + I .r 





/(*) 



c 2 



1 + x 2 



(A17) 



where C ~ 0.577 is the Euler constant. Note that if the probability distribution is symmetric 
asymptotically, i.e., c\ = C2, then /3 = 0, and a in Eqs. (|A15|) and ( |A17|) is the mean 
value calculated in the principal value sense. If the probability density is fully symmetric, 
f(—x) = f(x), then f(u) is real, a — (3 — 0, and the behaviour of the characteristic function 
at small u is especially simple: 



1 — cM 



After establishing the form of of f(ui) at small u, Eq. 
derive the limit theorem, starting from Eq. 



FJu) 



(A18) 

for a ^ 1, we can proceed to 



| c(l + z / 3sign(o;)tan^)n 1 - 



exp 



n 



-c ( 1 + if3 sign(u;) tan ^ J n 1 



(A19) 



for large n (this formula appears in the theorem by A. Ya. Khintchine and P. Levy as 
a canonical representation of stable probability distributions, see Ref. 0). Using the last 
expression in F n (X) = j- f e lujX F n (uj)(ko we obtain the limit distribution in the following 
form: 



F n {X) = n 



a c a fa/3 n a c a (X — a) 



where 



fat 



e iu,x-\u\* exp 



-ip sign(u;) tan — \oj\ 



2t7 



(A20) 



(A21) 



is a universal function of the two parameters, a and (3, normalized to unity: J f a (x)dx = 1. 
The results for a = 1 are obtained in a similar way, with a replaced by a + c-/3lnn in 
Eq. ( |A20| ) , and f\p(x) given by Eq. ( |A21| ), in which tan ^ is replaced with — - In \uj\. 

Equation (|A20| ) shows that for < a < 1 the limit distribution of the average widens 
with the increase of n, i.e., fluctuations of the average increase with the number of variables 

ot— 1 

averaged. Since n~^~a — > for n — > 00, the shift of the distribution ( |A20| ) by a is actually 
unimportant in this case and one can put a = 0. For a = 1 the shape of the distribution 
F n (X) does not depend on n, i.e., fluctuations are neither enhanced nor suppressed by 
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averaging. If (3 ^ the whole distribution is gradually shifted proportionally to In n into the 
direction determined by the sign of /3. For 1 < a < 2 the distribution of the average does 
become narrower with n, however the rate of suppression of fluctuations, X oc n~( a ~ l ^ a 
is slower than the standard CLT n~ x l 2 . Again, for symmetrically distributed Xi, the limit 
distribution F n (X) is even simpler, as a = (3 = in Eqs. ( |A20| ) and ( |A21| ). 

There are a few cases where f a p and, hence, F n (X), are known explicitly. For a = 1, 
(3 — (ci = C2) Eq. ( |A21| ) gives the Cauchy law, 



/i,o(^) 



1 1 



7T 1 + X 2 



(A22) 



and c = |(ci + c2). For a = 1/2, f3 = the limit function can be expressed in terms of the 



error function = 2tc x l 2 JJ* e * dt [IT]: 



|,o W 



20F|x| 3 /2 



Im 







1 


< e 4 e 4 * 








\2VixJ_ 





(A23) 



For the same a = 1/2 in the maximally asymmetric case, ci = 0, C2 > 0, i.e., /? = 1, which 
takes place if the random variables Xi are positive, one easily obtains the following simple 
answer IBJ: 



0, x < 0, 



(2n) 



-1/2x^-3/2 



x > 0, 



(A24) 
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